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A MONOCOQUE WING- MODEL LOADED IN BENDING- * 
By E. Schapitz, H. Feller, and H. Koller 



Bending tests with, transverse loads and with pure 
"bending were undertaken on a doubl e-web monoco que wing 
model in order to establish the relation "between the state 
of stress and the results from the elementary "bending the- 
ory. The longitudinal stresses in the stiffeners were 
measured with tensiometers and the shear stresses in the 
sheet were calculated from them. The measurements were 
made at "both moderate loads with no "buckles in the cover- 
ing and also at loads with which the critical stress in 
the individual panels was exceeded. For the comparison, 
the wing skin was considered as stiffened sheet according 
to the shear panel scheme. In this way the statically in- 
determinate calculation was confirmed by the test results. 



I. STATEMENT OF TEE PROBLEM 



The well-known type of aircraft construction in which 
not only the shearing forces but also the direct forces 
resulting from bending are resisted by the stiffened sheet 
covering is designated in a larger sense as the monocoque 
type. (See fig. 1, reference 1.) Such a construction of- 
fers advantages in respect to aerodynamic and weight con- 
siderations as compared to spar wings with torque tubes. 
TJie aerodynamic advantage consists in the fact that a wing 
shell stiffened by members in the di recti on of the span i s 
generally stiff er in regard to buckling than a torque tube, 
which is stiffened only by a single auxiliary rib. There- 
fore, in the case of monocoque wings there is less danger 
that permanent buckles might occur in the covering and that 
the drag might increase (reference 2). The weight savings 
are effective only if at least one side of the wing can be 
built without openings (for tanks, engine nacelles, and the 
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like) and if the questions of fastening io the fuselage and 
the connection "between the units, are solved satisfactorily. 
In this case it is especially important whether the wing 
construction in the neighborhood of the fuselage must be 
transformed into a spar construction. With shoulder cov- 
ers and deep covers, for example, this can be given up and 
one can attach the individual stiffeners at the connections 
by means of screws. If, hpwever, only a single strut is to 
be -carried through the fuselage (for example, at the prin- 
cipal span), then the uniformly distributed longitudinal 
force in the wing shell must be transmitted to this point. 
The same problem arises as soon as the shell is interrupted 
by an opening. 

For the investigation of the strength, one can resolve 
the questions arising in the. monocoque wing into a series 
of individual problems* The infliience of the single moment 
introduced at the point of attachment can be represented 
in such a. way that one superposes the states of stress ac- 
cording to the ordinary bending theory and a group of 
forces in equilibrium (fig. l). This equilibrium group de- 
creases with increasing distance from the point of attach- 
ment. Another individual problem is the increase of stress 
in the stiffeners of a wing loaded by a transverse force 
at an opening if the deformation of the same is not re- 
strained (fig. lc). In these cases the longitudinal stif- 
feners in the interior are not active in the neighborhood 
of the application of the load on account of the absence 
of shear forces at the free edge; while, according to the 
simple theory of bending of beams, they should show a higher 
stress- than the stiffeners at the spars because of their 
greater distance from the neutral axis. Figure lc shows a 
third problem; Because of the shear flow* the built-in 
transverse cross section would warp out of its plane if it 
were not restrained by the supporting forces. Through this 
resistance to deformation, there occurs a group of normal 
forces in equilibrium whose effect vanishes at the free 
end. Similar disturbances arising from the restrained 
warping of the cross section are present in the case where 
the load is distributed along the length, ■ 

The tests on a monocoque wing model described in this 
paper are a basi s for the experimental de t erminat i on of 
the stress disturbances. The variation of the actual 

*The term "shear flow" literally translated from "schub- 
fl&sse," appears, to be used rather extensively in 'the 
German literature to denote shear stress times sheet 
thickness. 
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stresses from those obtained "by the. elementary theory of 
"bending -are - taken as a measure of the stress disturbance. 
By introducing pure "bending moment from couples, the stress 
disturbances from the application of concentrated loads 
were studied. By introducing concentrated loads at webs 
and di st ribut ed loads over the length of the shell, the 
effect of the restrained deformation of the cross section 
was studied. The measurement s were therefore made at loads 
"before as well as after buckling. The shear stresses were 
determined from the test results in a manner similar to 
that carried out for a stiffened circular cylinder in an 
earlier investigation (reference 3). 

The individual problems studied can be treated theo- 
retically also with the methods for the calculation of 
■statically indeterminate systems by making simplifying as- 
sumptions. For the stiffened sheet, of which the monocoque 
wing is made, the process of calculation is carried out in 
a work by Ebner and Koller (reference 4). In it the sheet 
field is considered to resist shear only and the normal 
forces are considered concentrated in the longitudinal and 
transverse stiff eaers. The strength of the skin is taken 
into consideration by the addition of an effective width 
to the cross section of the stiff enor. In the case of 
tests of a circular cylinder loaded at a single point (see 
reference 3), it has been established that the calcula- 
tions, according to this "shear-field theory," give a very 
good agreement with the measured values for stiffened com- 
plete shells (reference 5). The stress measurements pre- 
sented herein can only serve to prove the appli cabili ty of 
the shear-field theory to the calciilation of stiffened 
sheet. In section IV the procedure of calculation is shown 
and a comparison between the calculated and measured stress- 
es is made. 

- II. TEST SPECIMEN, TEST PROCEDURE, AND LOADING CONDITIONS 



The test specimen, which is described briefly in an 
article by H. Ebner (see reference 1, pp. 185-186*) repre- 
sents the interior section of a double -web monocoque wing. 
Its statical tension side (under side) is, in the case of 
this first model, concentrated into two heavy chords on the 
webs ( f ig . 2 ) . The model rep res en ts therefore a t rough , 
whose open side is turned upward because of the greater 

The model and the arrangements for loading and supporting 
were built by Dr. Ing. Lucker. 
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convenience in the work. While testing, the model is 
therefore reversed from the position in the actual wing 
(compressive side down). The model, as an open trough, 
facilitated the stress measurements appreciably; it is 
used in similar form also in the case of American investi- 
gations (reference 6 ) . The tension chords , con si sting of 
formed duralumin sheet, each of 15 cu 2 (2.325 sq. in.) cross 
section, are attached with s crews to the angle s of the 
spars, in order that a tension skin could he used in its 
place lat er . 

The compression side consists of a sheet covering 1.2 
mm (0.0186 in.) in thickness, with a radius of curvature 
of 4»000 mm (157.5 in.). It is stiffened with closed hat- 
shaped stiffeners spaced 140 mm (5.51 in.). The clear 
width of sheet "between the rivet rows is 90 mm (3.54 in.). 
The center stiffener is composed of two C-shaped sections. 
At the edges, the covering is joined "by two heavy angles 
to the spar webs. Figure 3 shows the dimensions and cross- 
sectional areas of the^e longitudinal stiffeners. The av- 
erage sheet thickness (total cross section without spar 
angles, expressed in units of arc length) of the compres- 
sion side defined "by the inner rows of rivets at the spar 
angles 1 and 9, is 2.75 mm (0.082 in.), thus the' cross sec- 
tion is distributed 43.5 percent to. the skin and 56.5 per- 
cent to the stiffeners. 

The spar web is 3 mm (0.118 in.) thick and so propor- 
tioned that the spar can remain flat until failure in the 
case of an ultimate load test with transverse "bending 
loads. In the region of the fixed end, the web is provided 
with longitudinal stiffeners in order to increase the re- 
sistance to buckling (reference 7). 

The position of the centroid S of the whole cross 
section is indicated in figure 2. The moment of inertia 
of the cross section roferred to the transverse axis through 
S is J = 24,300 cm 4 (583. 8 in. 4 ). The distance z of 
the centroid of the stiffener from the transverse 

axi s through S (see fig. 13) is given in table I . 

As transverse stiffeners, there are cross frames 
spaced a distance of 500 mm (14.69 in.). The cross-frame 
chords consist of hat-shaped sections of 1.5 cm 2 (0.233 
sq. in. ) cro^s section; the vozlz and diagonals are C- 
shaped sections of 2.3 cm 2 (0., r 5o7 sq. in.) and 0 C 8 cm 2 
(0.124 sq.. in.), respectively. The transverse frames are 
attached by special angles at the spar stiffeners. At the 
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spars they are riveted so that their posts together with 
the wet) stiff ener lying on the outside, stiffen the weh 
sheet. The designation of the spaces and the transversa 
joints is shown in figure 2. 



TABLE I. Distances, z, of the Centroids of the 

Stiff en ers, Sg , from the Transverse Axi s 
through the Oentroid of the Section, S 



Stiffener i 


1 and S 


2 and 8 


3 and 7 


4 and 6 


5 


10 and 11 


a ) zi = 


14. 7 


15. 9 


17.1 


17.8 


17.7 


16.9 .cm 




5.79 


6.26 


6.73 


7. 01 


6.97 


6.65 in 




14. 8 


16 .3 


17.5 


18.2 


18.2 


16.9 cm 




5.83 


6.42 


6.89 


7.17 


7.17 


6.65 in 



a) S g refers to the stiffener with the strip of sheet 
"between the rivets. 



"b ^ Sg refers to. the stiffener with a strip of sheet 
equal to half the spacing on each side. 



For support ing the mo del, use is made of a frame (fig. 
which is self-contained and can support "bending loads 
without the aid of a foundation. The spar chords and the 
skin in the compression zone are attached through an an- 
gle at the edge to a heavy steel plate. Wedges are used 
to transfer the compressive force from the longitudinal 
stiffeners to the steel plate. By this means, a uniform 
"bearing of the stiffeners on the. plate can he obtained. 
(See fig. 5 for the attachment descrihed.) The compression 
sheet is reinforced directly at the fastening "by a fish- 
plate in order to exclude as much as possihle a failure di- 
rectly at the point of application of the load. The two 
tension chords are carried as fishplates through the heavy 
steel plate and fastened "by means of turned holts to a "box 
"beam arranged "between the supports of the loading frame. 

For the loading of the model, it is provided that at 
each weh "both a concentrated load at the free end and also 
a smaller transverse load can "be applied at each cross 
frame. Further, however, pure "bending moments in various 
planes can he introduced at the free end and, indeed, at 
"both wehs and at the center of the model. For applying 
the moments, use is made of a loading device (which has 
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"been described (reference 8) previously), which consists 
of a loading triangle and an attached rectangular link 
(fig. 4). The tensile forces are transmitted to the con- 
centrated tension chords by means of a transverse member 
(see figs. 6c and 6d) while the compressive forces act 
through a loading beam to the group of stiffeners. In the 
case of applying the moments over the webs, stiffeners 1 
and 9 (fig. 6) are loaded and when moment is applied at 
the center of the section, stiffeners 4, 5, and 6 are 
loaded. Figure 5 shows the assembly for applying load at 
this 'point . For applying the distributed load at the cross 
frames, the web is carried "below the edge angle and is 
equipped with rings; the load is provided by dead weights. 
The concentrated loads at the free ends of the webs are 
applied "by a screw jack which acts through a dynamometer 
onto the web (fig. 4). 

The arrangement of the measuring stations on the trans- 
verse sections of the wing model and the st at ions on the 
stiffeners is shown in figure 2. Huggenberger tensiome- 
ters with a multiplication ratio of 1:1,100 were used as 
measuring instruments. The original gage length of 20 mm 
(0.787 in.) was increased in consideration of the small 
strains to he measured to 100 mm (3.94 in.) by means of 
attachments. Since temperature changes, especially "by the 
sun 1 s rays, have a great influence on the measurements, 
the temperature was controlled continuously by a thermo- 
stat. The greatest permissible temperature change during a 
series of measurements amounted to ±0.5° C. 

The tests related to the loading conditions shown in 
figure 6. The first case concerns the application of two 
equal transverse loads at the free ends of the webs with 
free longitudinal distortion of the loaded cross sections. 
In the second case, transverse loads whose total magnitude 
equaled the load applied at the free ends in the first 
case, were applied at the intermediate transverse frames. 
With this loading, the influence of the restraint to dis- 
tortion at the fixed end, according to figure 1c, should 
be studied. Since beyond the first load point there is an 
unstressed field, ther.e is also a restraint to distortion 
of the cross section in the first panel (fig. 6b). The 
third case concerns the appli cat ion of a couple at each 
web with free distortion of the cross section at the load- 
ed end and uniform support at* the built-in end. . (See fig. 
6c.) In the last case (fig. 6d) a single moment is applied 
at the center of the wing with the same conditions of dis- 
tortion as in the third case* 



N.A. C . A . Technical Memorandum No . 915 



7 



III. TEST RESULTS AND THEIR EVALUATION 

!• The Stress Peaks in the Pour Loading Cases "before 
the Buckling of the Covering 

The results of the strain measurements at the individ- 
ual measuring stations were interpreted into stresses with 
an elastic modulus, E, equal to 740,000 kg per cm 2 
(10,525,000 lb. per sq.. in.), which is an average value 
from various test specimens. Under the assumption of lin- 
ear distribution of stress across the section of the stif- 
fener, the stress at the centroid (ov) of each longitu- 
dinal stiffener was determined from thq values of stress 
at the individual gage lines. These centroidal stresses 
are shown in figures 18, 21, and 23 for comparison with the 
values calculated "by the shear-field theory. In order to 
give a clear picture of the test results, smooth curves- 
were drawn through the measured values which are shown in 
perspective in the following figures. 

a ) Application of a single load at the fr ee ends of 
iM-ieos.- In figure 7 the centroidal stresses, 0" L , of 
the longitudinal stiff eners are plotted on the wing span 
(coordinate x) and the wing width (coordinate y). The 
distribution given "by the elementary "bending theory is 
shown "by the light lines; stresses increasing linearly 
with x and decreasing with y from the center o\it to 
the webs be cause of the curvature of the skin covering . 
The figure clearly shows two variations in the courses of 
the curves: At the free end (loaded end), the increment 
of stress in the x- direction in the interior stiffeners 
disappears while in the edge stiffeners at the web the in- 
crement is greater than that according to the elementary 
bending theory. At the fixed end, the stress peaks are 
over the webs and a decrease in stress toward the middle 
is to be noticed, so that the curves In the transverse 

direction are curved opposite to those from the elementary 
theory . 

The variation at the free end arises from the fact 
that there, because of a lack of edge stiffening, the shear 
stresses in the skin must vanish. Since an increase in 
stress in a stiffener can result only by the occurrence of 
shear forces in the web or in the covering, then no increase 
in stress at the free ends of the interior longitudinal 
stiffeners is possible. The middle of the wing, therefore, 
is only a little effective in the first panel. The varia- 
tion at the fixed end is explainable, on the other hand, 
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from a restraint to distortion of the cross section: he- 
cause of the transverse forces, the cross section would 
assume the distorted form indicated in the smaller sketch 
of figure lc; since pi vane sections must remain plane at 
the built-in end, then there arises a group of forces 
which are in equilibrium. They consist of additional com- 
pressive forces at the lower spar chords and ten si le 
stresses in the middle. Their effect decreases in the 
longitudinal direction so that at the middle of the model 
(about 1,500 mm (59.06 in.) from the point of application 
of the load) the stress distribution corresponds to the 
elementary bending theory. 

b ) Loading_of_the_^^ oyer the 

length of the_ wing.- Figure 3 shows the stress peaks cor- 
responding to the case in which the transverse load is ap- 
plied, not as a concentrated load at the free end but as 
a uniform load distributed over the length of the wing* 
The total transverse load is the same as in case a). The 
stress distribution according to the elementary theory is 
parabolic and the diagram is made up of straight lines be- 
tween the cross frames; over the width of the wing (y- 
direction) there is an increase from the edges toward the 
center as was expected previously. Also in this loading 
case there are variations at both ends; at the unloaded 
free end, tensile stresses arise in the spar chords on the 
compression side and in the center the compressive stress- 
es are higher than according to the elementary theory. At 
the built-in end, stress peaks at the spar chords and a 
decrease in the stress at the center appear as in case a). 

The variance at the built-in end results from the re- 
straint to distortion of the cross section as in case a), 
only the effect is relatively greater because the same 
shear force and same secondary stress arising therefrom 
stand in contrast to a smaller bending moment. The vari- 
ance at the unloaded free end is opposite in direction to 
the variance at the same place in case a) and is brought 
about by the fact that beyond the last point of applica- 
tion of the load, there is an unloaded wing member which 
restrains the distortion of the foremost cross section. 
The accompanying internal state of stress is the same ex- 
cept for magnitude and sign as at the built-in end; com- 
pressive stresses are produced in the middle of the wing 
and t ensile stresses at the webs. 



c ) Lo ading by concentrated couples „ of „pure__b en ding 
moment.- Figure 9 shows the stress peaks for the case of 



N.A.O.A. Technical Memorandum Ho, 915 9 



pure "bending if the moment is applied at the loaded end as 
-• — a couple concent rat ed at the webs. The elementary theory- 
would give a constant stress over the length of the member. 
. A variance enters in this case only at the points of appli- 
cation of the forces since no shear forces are transferred 
and therefore the cross sections are free of shearing 
stresses at a sufficient distance from the point of appli- 
cation of the load. In consequence of this, the cross 
sections, free to distort, remain plane. 

At the loaded end, the compressive stresses are con- 
centrated in the spars at the points of application of the 
forces. The longitudinal stiffeners in the middle are 
nearly free of stress and near the edge also no increase 
of stress is shown since here the shear stresses vanish. 
The equalization to the elementary distribution of stress 
occurs within four panels ("between the cross frames).* 

An accordingly similar process occurs if a couple is 
applied at the middle of the wing (fig. 10). The compres- 
sive forces were applied to the three middle ribs 4, 5, 
and 6 by the loading arrangement represented in figure 5. 
The initial stress in ribs 4 and 6 is somewhat smaller 
than that in number 5 because of the flexibility of the 
loading bar. The equalization to the elementary distribu- 
tion in the case of applying the load at such an -unusually 
weak spot occurs in about the same distance along the wing 
as in the previous case with the application of the moments 
at the webs. 

2. The Effect of Buckling of the Skin 
on the Curves of the Longitudinal Stresses 

Since the pattern of the stress distribution must be 
calculated for the buckled as well as the unbuckled skin, 
the previously described stress measurement s were extend- 
ed into the range of buckling. The upper limit of the 
load range was therefore determined from a consideration 
of the strength of the stiffeners and the tendency to 
buckle; permanent buckles in the skin were to be avoided 
as much as possible because of the succeeding measurements. 
For the several loading cases , the loads at which the buck- 
les were first noticed were exceeded by the following 
amount s : 

Degree of 
Load case excess 

Transverse force at the free ends of the webs 2.25 
Pure bending, moment at the ends of the webs 2,00 
jj^re bending , ^moment in the middle 1.90 

j|C " " • — — ■ — ' " " r ' ' ■ — — — • — — 

This loading case corresponds to that studied by Lovett and 
Ho dee. (See reference 6.) 
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In the case of the di st ribut ed. transverse load the 
"buckling region was not investigated* 

The effect of the buckles on the stress distribution 
consisted in a diminution both of the compressive stiff- 
ness (the effective width) and the shear stiffness of the 
skin. This effect will he shown "by means of the measure- 
ments for the several loading cases. Thus the stresses 
were reduced to a 'oasis of a unit load in order to he able 
to ma.'te a direct comparison, Further, for clear plotting, 
the stresses were multiplied by the cross section of the 
longitudinal stiff cners (stiffener plus the adjacent sheet 
between the rows of rivets); thereby the force P x was 
obtained without consideration of the adjoining sheet panel. 

When the wing support s a concent rat ed load at the 
free end, the skin in the neighborhood of the built-in end 
buckles because the moment increases over the length while 
at the loaded end the sheet remains unbuckled. Figure 11 
shows the longitudinal stresses corresponding to the unit 
load of 1,000 kg (2,204.62 lb.) at the free end, plotted 
as a group of curves on the length of the wing. One no- 
tices that the variation at the built-in end is somewhat 
more marked after buckling than before. This results, on 
the one hand, from the decreased shear stiffness of the 
skin which leads to greater distortion and to greater sec- 
ondary stress; on the other hand, the smaller effective 
width in the buckled state counteracts this effect and 
leads in the case of equal forces to greater stresses in 
the stiffeners. That all effects act together can be seen 
from the fact that the longitudinal stresses in the stif- 
feners in the central portion of the wing after buckling 
are in part greater and in part smaller than before buck- 
ling.- An approximate mathematical representation of this 
interdependence will be shown in section 17. 

In the case of the loading with pure bending by cou- 
ples in the planes of the webs, only the strips of cover- 
ing directly against the chord angles of the spars buckled 
in the limited load range. The effect on the vanishing 
point was unimportant and not to be recognized in the scat- 
ter of the data points. (See, for example, fig. 11.) 
With the application of the couples on the three interior 
stiff en ers (4, 5, 6, fig. 12), a definite increase in the 
stresses is apparent in the latter in comparison with the 
state before buckling ( reduced to a unit load of 1 , 000 kg 
on a lever arm of 89 cm). This effect is clearly trace- 
able to the decrease in the effective width in the cover- 
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ing sheet; the load curves for the remaining stiffeners 
remain unchanged since the buckling was limited to the 
highly stressed sheet in the region from stiffeners A to 6.... 
Also in figure 12, the distribution of the data points fol- 
low the faired curves. 

3 # Shear Plow 

a ) _Q_omputat ions based on the me asur ed longitudinal 
stres ses .- The measured longitudinal stresses can he used 
to determine the distribution of the shear stresses in the 
load range before buckling . One obtains , thus , an in- 
sight into the shear loading of the . skin of the wing shell 
in the neighborhood of the points of application of loads 
and fixation. Furthermore, this evaluation is important 
in making the comparison with the calculations according to 
the shear-field scheme. 

In this way , the stiffness under axial loads of the 
stiffeners and of the sheet covering are assumed concentrat- 
ed in the longitudinal and transverse stiffeners; there- 
fore, the shear stress within a panel of the skin is con- 
stant and can change abruptly only at the stiffeners and 
the cross frames. If the transverse frames are not joined 
directly to the covering, then one can consider either the 
stiffness under axial load (through the lateral bending 
stiffness of the longitudinal stiffeners) divided uniform- 
ly over the skin and thus obtain a system with closely 
spaced cross frames, or one concentrates the shear stiff- 
ness of the adjacent half panel with the stiffness under 
axial load of the cross frames into a single imaginary 
frame joined to the skin. In the first case, the shear 
flow t = s changes steadily in each longitudinal strip 
and only at the stiffeners do abrupt changes occur; in 
the second case, the shear flow likewise changes abruptly* 
at the separate cross frames and is constant within the 
panel. The compari son ■ with the shear flow obtained in the 
tests shows then to what extent these calculated average 
values of shear flow approximate the actual distribution. 

Since the 'longitudinal stresses in the skin were meas- 
ured only at the changes in the section, it is necessary 
to make an assumption of the distribution of the longitu- 
dinal stresses over the cross section of the sheet in or- 
der to make a complete computation of the shear flow. It 
is approximate to try a variation within a single panel 
that i s linear between the st resses mea sured at the edges 
of the panel. A comparison of the external bending moment 
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with the resisting moment calculated from the measured 
stresses shows that this assumption of the effectiveness 
of the full width of sheet "between two longitudinal stif- 
feners in the load range "before Truckling leads to small 
differences from the state of equilibrium. Table II shows 
the distribution of the resisting moment for section Tip 
(near the application of the forces) which lies in the re- 
gion of greatest variation of stress. The moments of the 
stresses in the skin were calculated according to the above 
assumption concerning the longitudinal stresses. .It shows 
a satisfactory agreement between the external and resist- 
ing moments. 

TABLE II . . Compari son of the External Bending Moment 
with the Resisting Moment from the Longitudinal Stresses 

Measured at Section TIP 



Loading 
case 


Moment 
of .the 
1 ongi- 
tudinal 
forces, 

Px 

mkg 


Moment 
of each 
stiff oner 

mkg 


Moment 
of the 
web 

mkg 


— T 

Moment 
of the 
skin 

' 

mkg 


Sum of 
internal 
moment s 

mkg 


External 
moment 

mkg 


Coupl e s 
at the 
webs 


2, 225 


21 


287 


176 


2,709 


2,714 


Couple 
at the 
middl e 


2,084 


25 


126 


411 


2,646 


2 ,714 



Moment in mkg times 86.81 equals moment in inch-pounds. 



For the determination of the shear flow from the cen~ 
troidal stresses, 0" L , the longitudinal forces, P x , in 

the stiffeners (stiffener plus the adjacent strip of sheet 
between the rivet rows) were determined and were joined by 
a faired curve over the length of the wing. The relation 

4J5s = Px . = t r - n ■ (i) 

obtained from the condition of equilibrium of an element 
of the member exists between the changes in the longitudi- 
nal force, P x , in a stiffener across the width of the 
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wing (coordinate x) and the shear flow, t^ and t r (con- 
cerning the designations , see fig ,13) .in the panels ad- 
jacent to the stiffener. By graphical differentiation of 
the faired P x curve, one obtains "by means of equation 
(l) the difference in the shear flow in the adjacent pan- 
els at each cross section and at each longitudinal stiffen- 
er. The plane of symmetry of the model through stiffener 
5 is also a plane of symmetry for the longitudinal stress- 
es in all the lo£>dixig cases. The shear stresses must there- 
fore always "be antisymmetric with respect to it. One ob- 
tains therefore "by means of equation (l) from one-half the 
T x curve for the middle stiffener, 5, the amount of the 
shear flow which occurs on each side of the longitudinal 
stiffener. If the four panels in one-half of the shell 
are designated 1,2,3, amd 4, the edges with the indices 
I and r in the direction of the stiffener numbering 
(fig. 13) and the shear flow at the lower edges of the webs 
with tg 



u 



then it follows from (l) that 



r X 5 


- 2 1 4r 






= t 3r - 




P 1 

r X3 


= 1 2r - 






= *ir - 


*al 




= *?u " 





(2) 



The distribution of shear flow in the sheet panels 
and with it the relation between t 4r and t 4 | ft t 3r and 

t 3 ^ , etc. , are found by the integration of the equation 

of equilibrixim of the elements of the skin 

— 5 + — = 0 
6 x 3y 

whereat a linear variation over the panel width, b, be- 
tween the edge values of a x ^ and cr xr is assumed for a x , 
therefore 

(y in each panel is measured from the left.) 
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to 



From this it follows "by differentiation with respect 



St 



_ _ ^x = . dcy xl ^ (j 

^ t riir ' \ 



dy dx v dx 

and. "by integration with respect to y 



dx 



T (x v) = T (x 0) - -^i y - (—21 12x1 ) Zf 
kX,y; T ^ X,U; dx 7 ^ dx ' dx / 2b 



(4) 



The shear stress is therefore parabolic over the cross sec- 
tion in each panel. The difference in the shear flow be- 
tween the edge.s I and r of a panel is 



t r - t l = 



IS (i£kE 

2 v 



d £si ) 

dx J 



ax 



(5) 



The terms 



dcr^ 



for the individual panels are 



and 

dx ax 

obtained from the measured stresses in the skin at the 
stiffeners which are joined by a faired curve and then the 
derivatives are determined graphically at the sections 
studied. Now if a summation is made, from the middle out 
to the. e;dge. at . a certain cross section, of the differences 
in shear flow calculated according to (2) and (5), a curve 
for the variation in the shear flow is obtained for the 
cross section under consideration. It consists of paraboli 
segments over the free width of the sheet according to equa 
tion (4); at each stiff ener there is a discontinuity accord 
ing to (2). If the skin stresses on both sides of the stif 
fener agree, then the adjacent- parabolic segments have the 
same slope. The discontinuity at the stiffener can then be 
divided into two separate discontinuities of equal -magni— ' 
tude at the rows of rivets,, whereat the shear flow in the 
strip of skin under the stiffener varies as a straight line 
in the direction of the ends of the tangents of the ad- 
jacent parabolic segments. (See figs. 14 to 17.) 



If 



'xZ 



is the force in the tension chord (fig. 2), 



then the shear flow, 
is, computed from 



at the upper edge of the web 



dP 



dx 



(2a) 



For the web (height, h, between the centroids of the 
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chords), equation (5) can "be rewritten in the following 
manner : ■ , . .. 



in which s< 



- t 



Su = 



h Sg | 1 dP xl 
2 L Pi d.x ■ 



1 dP, 



*z ax J 



(5a) 



is thickness of web 



3T X , cross-sectional area of the lower chord 
angle 

t cross-sectional area of the tension chord 



The relation (5a) 
tion of the shear 

and (4) 
lead to 



can serve as a check since the 
flow, tg^ , obtained "by means 

through the summation on the compressive 
the same result as obtained according to 



(5a) for the tension side. 



evalua- 
of (2) 

side must 
(2a) and 



"b ) The shear flow in the loading cases studied . - Th e 

shear flow in the loading cases studied was determined 

according to this process of calculation for the sections 

where stresses were measured. The results are shown in 

figures 14 to 16; the absolute amounts of the shear flow 

(for the measuring stations shown) are plotted over the 

half cross section. The check described for to is car- 

u u 

ried out in all cases; the maximum variation is 8 percent, 
which is traceable to the arbitrariness in the choice of 
the adjusted curve through the scattered data points and to 
the errors in the graphical differentiation. 

Of the four loading cases studied: bending with trans- 
verse forces (with concentrated loads at the free end and 
with loads distributed over the length) and pure bending 
(with forces applied at the webs and forces applied in the 
middle) - the last two cases differed from the first two 
in such a way that the elementary stress distribution gave 
no shear flow; therefore, shear stresses are traceable only 
to the distortions resulting from the application of con- 
centrated loads. The two cases of bending with transverse 
forces are different, in that in the case of the concen- 
trated load at the free end the total shear force to the 
absorbed over the width x remains constant while for the 
distributed load it increases toward the built-in end. 
The localized stresses at the points of fixation and at the 
application of the loads produce a distortion of the dis- 
tribution of the shear flow. 
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In the case of "bending with concentrated loads at the 
free end (fig. 14), the total shear flow in the web is 
concentrated, at first, "behind the point of application of 
the load "because the longitudinal stiffeners (except 1 and 
9) are not active. After that, the shear flow in field VI 
ri se s very rapi dly "because this field must t ransf er the 
total shear for the stress increments in the interior stif- 
feners. With the approach of the longitudinal stresses 
to the linear variation, the shear flow decreases over the 
whole width. The distortion at the fixed end causes a 
strong decrease in the shears over the whole width of the 
wing (see left side of fig. 14) "because the "buckling stress- 
es unload the inner panels and st if f eners so that the in- 
creases in the forces and the accompanying differences in 
shear flow from field to field are small er . 

In the case of bending with distributed load (fig. 15) 
the increasing transverse force corresponding to the shear 
increases from the free end to the fixed end until the dis- 
tortion resulting from the "buckling forces at the fixed end 
again cause a decrease which is more noticeable in the in- 
terior panels than in those near the webs. The increase in 
the normal stresses in these panels is (see fig. 8) notice- 
ably greater than in the interior panels. The reverse 
changes in the shear curves for cross sections la and 
13 at the "built-in end are caused "by the fact that the in- 
dividual "bending of the longitudinal stiffeners at these 
points unloaded the hacks of the stiffeners and the sheet 
stresses increased accordingly. 

In the case of pure "bending with the forces applied 
at the webs (fig. 16), the shear reduces to zero with ap- 
proach to the built-in end; it fades out steadily at all 
sections from the web toward the center. In the case of 
the application of the load at the center (fig. 17) , the 
shear attains its greatest value in the foremost cross sec- 
tion in the panels which are directly outside the loaded 
stiffeners, 4 and 6. At a greater distance from the points 
of application of the load, these maxima flatten out and 
the shear flow vanishes at the built-in end. 



N.A.C.A. Technical Memorandum No* 915 



17 



IT. CALCULATION OP THE STRESS DISTRIBUTION ACCORDING TO 
"' THE SHEAR FIELD ' CONCEPT " 

1* Basis of the Statically Indeterminate Calculation 

The stress distribution for the loading cases studied 
experimentally will "be determined mathematically with the 
assumption, that the skin resists only shear as in the pro- 
cedure developed in an earlier work "by Ebner and Koller 
(reference 4). This work is referred to later "by the ab- 
breviation E-K. 

The longitudinal forces in the stiffeners and the ad- 
jacent strips of sheet are considered concentrated at the 
centroids of the stiffeners;* the stiffness of the skin 
under axial loads is considered as an incree.se in the stiff- 
ness of the stiffener. Let it be assumed that the covering 
is not buckled at first. Then one-half the sheet panel on 
each side of the intermediate stiffeners is to be added. 
The web sheets are also regarded as being purely shear re- 
sistant with a height h - 30 cm (11,81 in.), and a thick- 
ness of sg = 0.3 cm (0.118 in.). In order to obtain the 
same stresses in the stiffeners at the webs of the ideal- 
ized system as in the actual model, it is necessary to add 
to that stiffener 'area an effective web area of j h sg = 
1.5 cin a (0.233 sq. in.). Under these assumptions, one ob- 
tains, with the dimensions from figures 2 and 3, the fol- 
lowing effective stiffener areas: 

F i = F 9 = 7.4 cm 2 (1.147 so. in.) 

F 2 = F 3 = F 4 = F e * F ? = F e = 3 - 8 cmG (°- 58 9 sq. in.) 

E 5 = 5.0 cm 2 (0.775 sq. in. ) 

F io = F n - 20 - 9 cm2 (3.240 sq. in.) 

The simplified system of stiffeners with these cross- 
sectional areas, which are joined only by a shear resist- 
ant sheet , occupies the same centroids as the actual model 
and their moments of inertia differ by only about 2 percent 
from the actual values. 

Let it be further assumed that the fixation of the 
model occurs at cross frame a ( s ee f i 2) since the 

The se longitudinal forces differ therefore from the earli- 
er values P x , which designated only the forces in the 
stiffeners directly in contact with the skin. 
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23.5 cm (9.25 in.) extension at the point of fixation is 
reinforced by the angle along the edge and a strip of sheet 
on the compression side. The end frame a is- assumed as ■ 
rigid. The transverse stiffness of the compression sheet 
i s composed of the stiffness of the skin to axial load in 
the transverse direction and that of the chords of the 
transverse frames, which are not joined directly to the 
skin - 1 1 i s assumed that the stiffness under axial load 
of the cross frame serves only to equalize the forces re- 
sulting from the curvature of the compression side and 
therefore the full transverse stiffness of the covering is 
effective in the shell. For the case of "bending with trans- 
verse forces, the normal forces in the transverse direc- 
tion and the transverse stiffness of the skin are consid- 
ered concentrated in the transverse stiff eners joined to 
the skin at the distance a = 50 cm (18.68 in.), with an 
effective cross-sectional area = 50 X 0*12 = 6,0 cm 2 

(0.930 sq. in.). In the case of pure bending, on the other 
hand, |g B = E 0.12 cm (E 0.0472 in.) is computed with 
the distributed transverse stiffness because for this prob- 
lem complete solutions of the equations of elasticity apply 
in the limiting case of heavy frames. 

The compression side of the wing model including the 
spar chords represents, under these assumptions, a stif- 
fened shell with nine stiffeners as they are treated in 
reference 4. The average spacing of stiffeners, b, is 
13.8 cm (5.43 in*); the skin thickness, s, 0.12 cm (0.0472 
in.); and the value of G/E is 0.385. Equation (18) of 
E-K then gives the following values: 



W L =1; 0,3 = = 298.7; «> S .= ^ « 3,760 (6 ) 



The orthogonal characteristic force groups of a 

regular shell with equal stiffeners are chosen as static 
redundant s; the partial forces are taken from table I of 
E-K for m = 9. They are given in table III, also the 
corresponding kuj values from which the shearing coeffi- 
cients of the elastic equations can be evaluated. The 
values of and Kg are taken from table I of E-K and 

multiplied by the corresponding values of cu from (6); the 
values of K L are determined exactly according to equation 
(19) of E-K. 
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because the attached stiff eners are not all alike. In the 
seating up of the elastic equations, therefore, "besides 
the coupling terms resulting from the strains in the. cross 
frames, there are also such terms, resulting from the lon- 
gitudinal stiff eners which are also neglected; that is, the 
values X of the redundant characteristic force. groups are 
determined from independent' 5-term equations, (See E-K , • 
example . ) On account of the curvature of the compression 
side, hearing stresses arise- at' the webs as a result of the 
characteristic force groups, as is hrought out in. section 
1-7 of E-K . Since the curvature is small, its effect in 
the loading cases with transverse "bending forces is neg- 
lected. With the introduction of couples , on the other 
hand, they must he considered because herewith the redun- 
dants assume rather large values. 



TABLE III 



Partial Forces and Shear Coefficients of the 
Statically Redundant Characteristic Force Groups 



Re- 
dun- 
dant 




Par t ial 


forces 


— 1 




Vis 




i =1 and 9 


2 and 8 


3 and 7 J 


^4 and 6^ 


5 


*ik 


-0. 940 


-0.5 00 


0.174 


0. 766 


1 


186 . 5 


2873 


~1936 


-2k 


. 766 


-.500 


-.940 


.174 




202 . 3 


813 


109 .6 


-3k 


-.500 


1 . 000 


-.500 


-.500 


1 


220.1 


448 


17.52 


~4k 


.174 


-.500 


. 766 


-.940 


1 


232. 0 


346 


4.34 



With the carrying out of the statically indeterminate 
computations, a few designations are changed from those 
used in the previous sections and also in comparison with 
the work of E-K. The transverse frames and the correspond- 
ing redundants are designated with digits 0, 1, . 6 
from the free end instead of with letters g', f, a; 
the index k indicates, generally, one of these frames. 
(See fig. 2, numbers in parentheses.) The stiff ener num- 
bering begins at the spars with 1, 2, the sheet panel 
between the stiffeners carries the number of the stiffener 
on the left (fig. 13). 
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2. Stress Distribution in the Case of 
Transverse Bending Forces 

Let the transverse stiffness of the compression side 
"be concentrated in the individual frames 0» l t ,.,,'"6, 
which are joined to the skin. Frames 0 to 5 have the 



same cross .section 



ss 6.0 cm 2 (0.930 sq. in.). (Frame 
6 ds rigid*.'} Since the extension at the free end, with a 
length of 23.5 cm (9.25 in.), is not closed "by a frame, 
then on the "basis of the shear-field notion- there can he 
no longitudinal stress in the field between the stiffen- 
ers and with it no redundant at the frame 0.- At frames 



1 to 6, characteristic force groups X 



produce the nec- 



essary internal state of stress in the stiffened shell, 



whose multiples X 



are calculated from the independent 

equations of elasticity. The coefficients of the equa- 
tions of elasticity are obtained from the values of ko> 
in table III, according to equation (17) of E-K , whereby 
it is considered that the rigid frame 5 produces no con- 
tribution, and the redundant X^ effective -there applies 

to the shear field only. It is then 



aE8 kk 


= 8k l Wl 


+ 




+ 


3 
2 




for 


k 


= 1, 2, 


3, 4 


aS6 e5 


= 8^0^ 


+ 




+ 


'5 
4 












aE8 66 


= 4^ 


+ 




+ 














aE6 k,k-i 


= SKl^ 












for 


k 


= 2, 3, 


4, 5 


aESes 






K 3 U>b 




1 

2 












aE5 k,k~2 










1 

4 


KgUJg 


for 


k 


= 3, 4, 


5, .6 



The 6 values multiplied by c = aE 10 
in tables IV to VI I . 



~ 2 



are assembled 



As a principal substitute system for the external 
loading , the 11 i deal beam" whose stat e of stress is ob- 
tained from the elementary theory of bending, , is taken as 
a basis- The internal state of stress of the characteris- 
tic force groups then giver the deviations from this sim- 
ple distribution which arise from the restrained distor- 
tion of the cross section at the built-in end and at the 
points of appli cat ion of the loads • 
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Loading case a , .cgn.ce.nt rat ed load at the free end' 
The . transverse load, Q, half to each web, acts at the 
free end. The bending moment at the ~k*k - frame amounts, to 

M k = (lea + 23.5 cm) Q 

M k = (ka + 9.25" in. ) Tl ■ 

With the distances to the, centers of -gravity of the stiffen- 
ers, zj[, in table lb, one obtains the distribution of 
force according to the elementary theory 

P <£> = - (k a + 23.5 cm) lA^A ft 

P ik } = - (k a * 9 - 25 ~j~ ^ 

and the corresponding shear flow, constant over the length, 

t (o) = - t<°> = - i - 5 > q 

4k 5k 2 J 

t (o) = - t (o) = . *_A!L_!_!i!i Q 

3k 6 k J 



et c . 

For Q = 1 kg (2.205 lb.), the values of shear flow 
multiplied "by J are given in table VIII. Since these 
shear flows are constant in the longitudinal direction, they 
load only the two end frames, 0 and 6. Table VIII gives the 
values of the forces Q> o ) multiplied "by j/b for frame 0 

at the ± ba atiffener resulting from Q = 1 kg (2.205 lb.), 
which are obtained by summing up the values J tl?' from the 

1 -lC 

edges. The corresponding values resulting from the redun- 
dants a-re also summarized in table VIII. (See section III 
of E-K. ) 



In this principal substitute system, the longitudinal 

forces make no contribution to the loading terms in 

the equations of elasticity because they are determined 
under the. assumption that .plane sections remain plane, so 
that through them there arises no distortion of the cross 
sections and therewith no additional work by the forces of 
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the static redundant s . * From the shear flow which is con- 
stant over the length there enters a loading factor differ- 
ent from zero only for the redundant X^ 6 at the built-in 
frame : 



aE =i Js G/E ^ ik 



Since frame 6 is rigid, then only the end frame (0) yields 
loading factors and, indeed, only for the redundant X 

-(H) 2 t> 3 I J 0 (o) / AP (H) . R (tO\ 




If one inserts the numerical values of the dimensions as 
given previously and uses table VIII, then the loading 
terms (multiplied by C = a E 10~ 2 ) contained in tables 
IV to VII inclusive, are obtained for Q, = 1 kg (2.205 lb.). 

The solution of the 5-term equations is carried out ■ 
according to a usable system and the results are given in 
table IX, One recognizes from the equations of elasticity 
that the transverse stiffness is of importance only for the 
characteristic force group X jj. , and further, that these 

redundants are appreciably greater than the other groups. 
In this way, the disregarded coupling with the force groups 
of higher order is small and is nonessential for the final 
force di st ribut i on . 

By superposition of the values for the redundants from 
t abl e_I X , one obtains the force 



^Because of the curvature of the sheet, the summations 

2 — were not equal to zero. If one considers, 

i = i E -B ± 

however, the effect of the supporting 'for ces , the charac- 
teristic force groups in the curved sheet occur in the at- 
tached spars (see section IV of E-K) , so the contribution 
of this work vanishes. 
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P P (o) + 2 P (p,) X 

and, by division by the effective cross section of the 
stiffener, one obtains the centroidal stresses in the stif- 
feners at the cross frames, .In figure 18, these stresses 
are plotted along the length of each stiffener for the 
load Q, = 1,700 kg (3,749 lb.). According to the idealized 
model with the cross frames attached to the skin, they must 
be joined "by straight lines. Since the actual frames are 
not 'located in the skin, then the shear in each longitu- 
dinal strip changes regularly and the values of the stress- 
es at the frames are therefore to "be adjusted "by a smooth 
curve. For comparison, the centroidal stresses determined 
from the measurement s are plo 1 1 ed in figure 18 . It can "be 
seen that the di st o rt i on of the elementary stress distri- 
bution is well represented "by the statically indeterminate 
calculations . 

The distribution of shear flow follows "by superposi- 
tion from table VIII , and the values for the redundants: 

*ik = t ik^ + j Jfc | 1 t ik ) (X p,k ~ X n f k-i* 

The stepped curves for the shear flow in the individual 
longitudinal strips calculated in this manner for Q, = 
1,700 kg (3,749 lb.), are adjusted by the continuous 
curves with the same integral value and therefrom the val- 
ues of shear flow are obtained for the various sections 
where measurements were made. They are shown in figure 14. 
Because the longitudinal stiffness of the skin is concen- 
trated in the stiff eners, the shear flow in the transverse 
direction remains constant within a panel. By comparison 
with the shear flows which were obtained from the measure- 
ments, one sees that these constant shear flows, in gen- 
eral, show the same distribution. Since the determination 
of the shear flow from the measured value sis subj ect to 
inaccuracies because of the scatter, so the deviations 
are much great er than in the case of the compari son of the 
measured and § calculated longitudinal stresses. 

From the preceding calculations, an estimate can be 
made of the effect of the buckling of the skin at the 
fixed end on the stress distribution. The loading terms 
at the fixed end increase with a dec r ease in . the shear ■ 
stiffness sG ; the coefficients in the equation, however, 
do not increase in the same rate because in the contribu- 
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tion from the stiffeners and .frames the smaller effective 
widths of skin cannot work out so rapidly. Consequently, 
the values of the redundant s become larger - that is, the 
distortion of the elementary distribution must "be more 
marked at the fixed end after the buckling of the skin, as 
is confirmed in figure 11 "by the measurements. 

Loading case ^ distributed transverse load on the 

webs : At the transverse frames 1, 2, 3, 4, and 5, a trans- 
verse load q is appli ed , hr.lf to each spar. The bending 
moment at the k"^ 1 frame amount s to 




The shear flow in the k uh transverse strip, accord- 
ing to the elementary bending theory, for* q = 1 kg 
amounts to 

( 0 ) 1 ^5*5 

4k 2 J V a J 



( 0 ) = ^ \ Ps^b tJ± Z ± (3S„Z 

3k ~ " " J K a 

et c . 

These differ from the former values in loading case a) 
only by the factor l/a (M^ - M^- x ) . For the redundant s 
at each intermediate frame, loading terms arise from the 
shear in the sheet and frames, which can be combined with 
the former value s on the basi s of relations between the 
shear flows. If we designate the former values with the 
index E, then we obtain (for q = 1 kg) a contribution 
from the shear in the sheet. 

. ( 6so)e (i'trJi- . ik ±J _:JV) 

= - a (K k-i " 2a k + M k+i> <'s«'l 
and after introducing the above expressions for M^: 
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for - k 



2, 



4, 5 




+ 5 ( S 60-)s 



For the contribution from the Gross frames, one finds, 
correspondingly : 




= 7 (M. - 4M V , + 61L - 4M. ' . + M. ) (S 10 )™ 
a k-2 k- 1 k k+ 1 k+s 1 ° 



This value differs for k = 2, 3, 4;'" in comparison, it is 



The loading terms, so determined and the corresponding so- 
lutions of the e quat i ons of elasticity' for q =.1 kg are 
contained in tables IV to VIX and IX. 3Por the load q = 
340 kg (750 lb.),' the final stiff ener stresses calculated 
by superposition are plotted in figure 18 and are compared 
with the average" values from the measured stresses. The 
distribution of shear flow, determined as in the loading 
case a) is shown in figure 15. Also in this loading 
case, the results of the statical indeterminate calcula- 
tions are satisfied by the measurement s • The variations 
at the built-in end rest partly on the variations in the 
measured stresses and partly on the neglected coupling of 
t'he redundants in the statically indeterminate calculations. 



The departures from the elementary stress distribu- 
tion in these loading cases (figs. 6c and 6d.) arise through 
the introduction of concentrated forces in the 'stiffened 
sheet of the compression side. At the point of application 
of the lo.ad, the departure is recognized : as the difference 
between the actual and the . elementary distribution of forces. 
This equilibrium group can be received into groups of the 
form of the orthogonal redundant s whi ch have , as a re suit , 
independent of one another, a decreasing curve of the cor- 
responding redundant in the interior of the shell. The de- 
crease can be calculated for each redundant according to 
equations solved in section 111,3 of E-K. Therein, if the 
transverse stiffness is assumed uniformly distributed (lim- 




3. Stress Di st ributi on in the Case of 
Pure Bending Moment Introduced -by Goupl es 
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itin? case of cross frames spaced infinitesimal disto-nces 
apart), one obtains a continuous curve of characteristic 
force groups dependent on the longitudinal coordinate x 
(calculated, out from the loaded cross section), 

Next, the curve of the redundants is calculated if the 

equilibrium groups Fx , T 2 , T 3 , F 4 from the components Q f 

the characteristic force groups X , X 2^ 3 ^> —4k ( see 

table III) act at the free end (x = 0). For the charac- 
teristic force group X x 

*- L = ~ = 0.2238; K s o) s = a 2 K s w> s = 4.84 X 10 6 

therefore, according to equations (28 l ) and (29') of E-K : 

1 = / ± a (0.0320 - = / e 0.01766 

p y * s w s y ^s^s ^0.0128 y 

With this one obtains, according to equations (27 1 ) and (29 ! ) 
of E-K, the curve of the redundants for the application of 
the equilibrium group I\ at the end cross section: 

X, (x)=e~ °- 0320 X (cosh 0.0128 x + 2.50 sinh 0.0128 x) (7) 
in the case of frames elastic in the transverse direction, 

X x (x)=e~ °- 0320 x (cosh 0.0128 x + 1.45 sinh 0.0128 x) (7«) 

in the case of a rigid frame at 0 and elastic intermediate 
frame s , 

X,(x)=e- 0.01766 x (7 „) 

in the case of rigid frames. The curve of the characteris- 
tic force group of higher order can be computed under the 
assumption of rigid frames. (See E-K, example no. l). 
Thus one obtains : 
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o . 034 6 x 



X 2 (x) =r e 
X 3 (x) = e 
X 4 (x) = e 



O . 0567 X 



O . 04 8 6 X 



• ••■ ... . . ( 8 ) 



The curves for this characteristic forcie group according 
to these formulas are plotted against the shell length in 
figure 19. The values "become negligibly small at the 
built-in end of the shell (x = 323.5 cm (127.3 in.), so 
that the above equations valid for infinitely long shells 
can he used. 

XfQg;M^g-Case„c, L _ couple s, applied to the spars : Pi gur 
20 shows the forces introduced in the chords of the spars 
for a load less than that producing "buckling. The total 
"bending moment introduced amounts to M B = 2,714 mkg 
(235,500 in. -lb.). The forces P{ 0 ' of the elementary 
distribution are obtained from 



with the values of z i from table I,b and J = 24,300 cm 4 
(584 in. ). The equilibrium group, which represents the 
difference between the two distributions, is* resolved into 
an t equilibrium group r, whose longitudinal forces in the 
stiffeners of the compression side form within itself an 
equilibrium system, and small additional moments in the 
chords of the spars which are in equilibrium with the su- 
perimposed forces of the curved compression side. (See 
fig. 13 of E-K . ■) 

r (with the partial forces P^°^) are resolved into 

the equilibrium groups IV to T 4 from the partial forces 
( M* ) 

P i of the characteristic force groups . X„: 



The development coefficient c is obtained regarding the 
orthogonality relation ^ 




r = 



c i r i + c 2 r a + c 3 r 3 + c 4 r 4 



(9) 



9 



(M-) p (v) _ 




0 
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for \i ^= v according to. eauation (26) of E-K 

, -< r 

(.1 *M ,[» . |) ' 

After inserting the numerical values from figure 20 and 
table III, it is found that 

c x = +■ 1,528; c 2 = - 1,103; c 3 = + 800 5 c 4 = - 207 

Let it he assumed that T exerts the principal ef- 
fects of the equilibrium group; the additional moments 
from the curvature of the shell fade away over the length 
of the shell approximately as the characteristic force 
group of the first order. The • redundant s resulting from 
r are obtained by multiplication of the values in figure 
19 with the corresponding constant c^. For X l5 the val- 
ues from (?) (elastic transverse frames) are to be used 
because at the free edge of the loaded side no transverse 
stiffness of the compression skin is present. If one su- 
perposes now the elementary distribution of the longitudi- 
nal forces resulting from these redundantts, then the final 
distribution of forces is obtained, and from this the cor- 
responding stress relations which are shown in figure 21. 

The comparison with the test results gives a good 
agreement for stiffeners 3 to ?. In stiffeners 2 and 8 
the measured stress in the first panel is not as high as 
obtained by computation. The stresses in the chords of 
the spars fade away correspondingly more slowly. These 
variations rest for the most part on the coupling of the 
redundants by the nonuniform cross sections of the stif- 
feners, because the characteristic force groups used are 
orthogonal only with equal stiffeners. If one sets up the 
orthogonal characteristic force groups belonging to the 
actual cross sections of the stiffeners of the model ac- 
cording to the general procedure outlined in section 
111,2 of E-K, and carries through the computations, then 
one obtains the stress curves shown as dot-dash lines in 
figure 21, which show good agreement with the measured 
values • 



For the calculation of the shear flow at the. cross 
sections where measurements were made, one first determines 
the values 
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dx 

at these places whereat one ohtains the derivative 'In 

"by differentiation of (V) and (8)' and multiplication with 
the corresponding constant c^. Then 

. , . 1 d? 4 . 

t3(x ) = + 

* dx dx 
etc. 

The results of the numerical calculations are plot- 
ted in figure 16. One obtains essentially the same dis- 
tribution of shear as given by the measured values, to "be 
sure, with greater variation's as in the loading cases with 
transverse "bending forces. 

Loading case &, coupl es allied at the middle of the 
2A£5 • The couple had the same moment M = 2,714 mkg 
(235,500 in . -lb • ) as in the previous loading case; there- 
fore, it gives the same elementary distribution of forces. 
The compressive force of 7,050 kg (15,545 lb.) is applied 
through the attachment in figure 5 to stiff eners 4, 5, and 
6, whereat the distribution to the individual stiff eners 
depends on the flexibility of the attachment. In order to 
avoid the determination of this distribution from a stat- 
ically indet erminant calculation, it was assumed as corre- 
sponding to the results of the measured stresses in cross 
section VII a 1 . (See fig. 2.) One therefore obtains the 
forces given in figure 22 eit stiff eners 4, 5, and 6. 

The method of calculation' is exactly the same as for 
.the . preceding loading case. As "development coefficients 
of the equilibrium group r, one obtains the values : 

c x = - 1 ,296 ; c 2 = - 700; c 3 = 157; c 4 = + 301 

For the fading out of. X x , equation (7 1 ) (rigid end frame) 
is chosen because an appreciable stiffening' in the trans- 
verse direction is produced by the attachment at the cen- 
ter. The additional moment' in.' the spars resulting from 
the curvature of the sheet should decrease* as the force 
group X x in the longitudinal direction. 
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Figure 23 contains the result s of the numeri cal calcu- 
lation for the stresses in the indi vi dual longitudinal 
stiffeners. Also here there is good agreement with the 
measured values. The corresponding distribution of shear 
flow is shown In figure 17; the stepped curve of the calcu- 
lated shear flow over the cross section furnishes thereby 
useful average values for the shear loading of the individ- 
ual stiffeners, 

V. SUMMARY 



By means of strain measurements in a hollow wing model, 
the longitudinal stresses in the compression side resulting 
from loadings with transverse loads and pure "bending were 
determined and plotted in perspective in the form of stress 
peaks on the outline. It was found that the calculations 
of such a framework according the elementary "bending theo- 
ry which neglects the effects of shear strains, may result 
in only a first app roximation. 

Through the shear deformations of the effective width 
of skin covering, distortions of the cross section arise 
with "bending produced "by transverse loads whose restraint 
at the "built-in end had as a result an unloading of the in- 
terior longitudinal stiffeners and a correspondingly higher 
loading of the chords of the spars , In the practically im- 
portant cases of distributed load over the length of the 
member, these superposed stresses can amount to more than 
30 percent of the elementary stress distribution. Small 
differences occur also at the points of application of the 
transverse loads . 

In the case of loading through pure bending by means 
of couples, concentrated forces. were introduced into the 
shell construction as they occur practically at the points 
of attachment of the shell wings. Also in this case,' the 
elementary stress distribution was distorted by secondary 
stresses. In the model studied, these distortions of stress 
extended in the longitudinal direction from the point of 
application of the load for a distance of about 1.5 times 
the width of the shell. 

Exceeding the buckling limit of the skin covering 
changed the stress distribution in the longitudinal stif- 
feners only slightly. In the case of bending under trans- 
verse loads, the distortion of the stress at the built-in 
end increased somewhat as a result of the decreased shear 
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stiffness of the skin. In general, the stresses in the 
longitudinal stiff eners increased "because of the decreased 
effectiveness of the skin with respect to longitudinal 
forces. 

From the longitudinal stresses measttred • bef o r e buck- 
ling, the distribution of the shear flow was calculated 
in order to obtain a picture of the shear loading in the 
skin introduced by that kind of stress distortion. 

The loading cases studied were treated theoretically; 
also , under the assumption that the skin "between longitu- 
dinal and transverse stiffeners resisted only shear. The 
carrying out of the statically indeterminant calculations 
was done according to the procedure established by Ebner 
and Koller in an earlier work ( Luf tf ahrt f or schung , vol. 15, 
1938, pp. 527-542). The values of stress determined math- 
ematically have good agreement with the measured values so 
that the "shear field theory" represents a usable idealiz- 
ation for the calculation of that kind of stiffened shell 
const ruction. 



Translation by Marshall Holt, 
Aluminum Research Laboratories, 
Aluminum Company of America. 
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TABLE IV 

Equations of Elasticity for the Redundants Xik 



k 


X l.l 


X a. . s 


X 1.3 


X !.4 


X !.6 


X 1.6 


Loading 
case a) 

<> 

0 


Loading 
case d ) 
nR ( l) 


1 


10.14 


-4.44 


0.48 








-0.155 


0.963 


2 


-4.44 


10 .14 


-4.44 


0.48 






0 


. 808 


3 


.48 


-4.44 


10.14 


-4.44 


0.48 




0 


. 808 


4 




.48 


-4.44 


10.14 


-4.44 


0.48 


0 


. 808 


5 






.48 


-4.44 


9 . 66 


-3.48 


0 


. 963 


6 








.48 


-3 .48 


4.11 


- . 808 


-4.195 



TABLE V 



Equations of Elasticity for the Redundant s Xgjj. 



k 


X 


2 . 1 


X 


2 . 2 


X 2 , 3 


S 2 . 4 


X 


2 . 5 


X 2 . 6 


Loading 
case a) 
( 2 ) 


Loading 
case o ) 
( 2 ) 
05 ko 


1 


3 


.409 


-0 


.518 


0.027 










0.0272 


-0. 1924 


2 




.518 


3 


.409 


-.518 


0.027 








0 


-.1652 


3 




. 027 




.518 


3.409 


-.518 


0 


. 027 




0 


-.1652 


4 








.027 


.518 


3.409 




.518 


0.027 


0 


-.1652 


5 










. 027 


-.518 


3 


.382 


-.464 


0 


-.1924 


6 












.027 


1 


.464 


1.650 

« 


..1652 


.8532 
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TABLE 71 

Equations of Elasticity for the Redundant s JCj^ - 



k 


X 3.1 


X 3 .3 


X 3 . 3 


X 3 . 4 


X 3 . 5 


X 3 .6 


case a) 


X . r\ a A ^ n £? 
Ju Q cL U.X il^> 

case h) 
• ( 3 ) 


1 


2.683 


-0.025 


0.004 








-0.0093 


0.0753 


2 


-.025 


2 .683 


-.025 


0.004 






0 


.0660 


3 


.004 


-.025 


2.683 


-.025 


0.004 




0 


.0660 


4 


— • 


.004 


-.025 


2.683 


-.025 


0.004 


0 


.0660 


5 






. 004 


-.025 


2 .679 


-.017 


0 


.0753 


6 








.004 


(-.017 

1 . i , 


1.333 


-.0660 


| -.3393 



TABLE YII 

Equations of Elasticity for the Hedundants J£ 4k 



k 


X 4.X 


X 4. 2 


X 4. S 


X 4.4 


X 4 . 5 


X 4 . 6 


Loading 
case a) 


Loading 
case h) 

06 k o 


1 


2.555 


0.114 


0.001 






1 


0.0025 


-0.0159 


2 


.114 


2.555 


.114 


0.001 






0 


-.0134 


3 


.001 


.114 


2.555 


.114 


0.001 




0 


-.0134 


4 




.001 


.114 


2.554 


.114 


0.001 


0 


-.0134 


5 






.001 


. 114 


2.554 


.116 


0 


-.0159 


6 






mm 


.001 


.116 


1.275 


-.0134 


.0695 
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TABLE VIII 

Shear Flow and Force in Transverse Frames Resulting from . 
Bending Under Unit Transverse Load Q = 1 k? (2.2 lh.) 
at Free End (load case a)) and the Static Redundants 







Shear flow 




Forces 


in transverse 


frames 




i 


Jt (0) 


at ik 


at (2) 


at (3) 


at ik 


b ^0,1+1 


£&,i+i 


R (2) 


( 3 ) 
%c,i+i 


P (4) 


1 


-243,1 


0.940 


-0.765 


0.500 


-0.174 


243.1 


0.940 


-0.766 


0.500 


-0.174 


2 


-181.2 


1.440 


-.2 66 


-.500 


.326 


424.3 


2.380 


-1.032 


0 


.152 


3 


-114.7 


1.266 


.674 


0 


-.440 


539.0 


3,646 


-.358 


0 


-.288 


4 


-45.5 


.500 


.500 


.500 


.500 


584.5 


4.146 


.142 


.500 


.212 



T&BLE IX 

S tat i c Redundant s Resulting from Bending 
"by Transverse Forces 



k 


Loading case a) 
Transverse load at the 
_ end of the wet) 
(Q = 1 kg (2.2 It). 


free 

) 


Loading case "b ) 
Distributed transverse load 

(q = 1 k? (2.2 lb. ) ) 




X lk 


X 2k 


X 7k 


X 4k 


X lk 


X 2k 


X 3k 


X 4k 


1 


0.0217 


-0.0082 


0.0035 


-0.0010 


-0.199 


0.065 


-0.028 


0.006 


2 


.0172 


-.0013 


.0000 


-.0000 


-.264 


.068 


-.025 


.005 


3 


.0233 


-.0002 


.0000 


.0000 


-.241 


.068 


-.025 


.005 


4 


.0526 


-.0014 


.0000 


-.0000 


-.083 


.062 


-.025 


.005 


5 


.1315 


-.0145 


.0003 


.0005 


.343 


-.005 


-.027 


.006 


6 


.3012 


-.1041 


.0495 


-.0105 


1.301 


-.519 


.254 

1 


-.054 
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&-3*i Long! tudinal stiff ener 
Ribs 

Test stations 



-Restraint » w « « ft w frLoad 

Test sections station 

Figure 2.- Transverse section and 
—3 longitudinal view of 

_j the wing model studied. Position of 



Distortion measurement stations. 
Sections 

Section 5 

J 



Secondary stresses 




U-it J 




at built-in section § 

(a) Sketch of monocoque wing. |" 

(b) Distribution of stresses in 
the case of loading 

by couples at the spar chords. H5Sm * 

(c) Loading by transverse forces figure 3.- Dimensions of the longi- 
at the spars. tudinal stiff eners on the 

Figure 1. Statement of the problem, compression skin. 

Wedges 




Longitudinal stiff ener 



Loading beam 




Plate 



Figure 5.- Loading of the compres- 
sion skin at the center. 



Loading beam 



H.A.C.A. Technical Memorandum Ho. 915 
Restraint _L**d- 



4.6,7.8 



Figure 4.- Test set tip. 




ft! JSiJJ** ^ co ** en \ r * ted transverse forces at the ends of the web. 
b Bending by uniformly distributed transverse forces on the webs 

(c) Pure bending by couples at the spar chords 

(d) Pure bending tqr a couple at the center of the wing 

Figure 6.- Loading cases studied. 

Stresses according to^v.^ Stresses according 

elementary bending 
theory 



Measured 
stresses 




Figure 7— Longitudinal stresses 

in the case of bending 
with a concentrated transverse 
force at the ends of the webs 
(Figure 6a v total transverse load 
force Q-1700 kg [3760 lb.]). 



Figure 8.- Longitudinal stresses 

in the case of bending 
by uniformly distributed loads on 
the webs (Figure 6b, total trans- 
verse load $ s 1700 kg [3760 lb.] ). 
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*J Stresses according to 

K elementary bending theory 

Measured 



f 




Jigs* 9,10,11,12,13 

Stresses according to 
plementary bending theory 



Figure 9.- Longitudinal stresses in 
the case of pure bending 
by couples at the spar chords 
(Figure 6c, total moment Mb « 2714 
mkg [335,500 in—lb.]), 

a a J 0 b off c a M fi d a^jf e a r /3 f (Z S j9^f' 




5? 5 




Figure 10. - Longitudinal stresses in 
the case of pure bending 
by a couple at the middle of the 
wing (Figure 6d t moment Mb -2714 mkg 
[235,500 in.-lb.j). 



« Load stage Mb = 2714 mkg 
(preceding buckling) 
— Load stage Mg —5160 mkg 
(buckling tone) 
Load stage 1700 kg 
(preceding buckling) 
- Load stage $ = 4500 kg 
(buckling sone). 

Figure 12.- Effect of buckling of the skin on 

the stress distribution in the case 
of bending byfa couple at the middle of the wing. 




Figure 11.- Effect of buckling 
of the skin on the 
stress distribution in the 
case of bending transverse 
loads at the free ends of the 
webs. 



Flange 



T 



Panel 



t 

L 



Figure 13.- Designation of the 
panel* in the. akin. 



V.A.C.A. Technical Memorandum Vo. 915 



Tig*. 14,15 



f 2 J i 

JL JtL Jtl Jl 


/ 4 

1 ± Jl 


i jf 


L JtL i 




















i 










"Hi 
























1 k 
















< 









^ — 1 


i 














t — i 


3 


1 1 








+ 




w 






i — i 


1 

rz 




































T— ■ 


a 1 A 1 
a \ \P 

' 1 1 ' 


V *f 

-r — H 


jr t 


eft- 








r^— i 1 1 

I — ! — i — I 


— # 


: :- 1 


i 
i 






= 


i i i , 


JM 




i 





Figure 14.-* Shear flow in the case of bending py concentrated loads at 
free ends of the webs (total transverse load, Q = 1700 kg 
[3760 lb.] heavy curve: shear flow determined from the strain measure- 
ments. Light horisontal lines: calculated values.) 




Figure 15.- Shear flow in the case of bending by uniformly distributed 
transverse loads on the webs. Total transverse load Q= 1700 kg [3760 
lb.] heavy curve: shear flow determined from the strain measurements. 
Light horisontal lines: calculated values. 
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rig*. 16,18 

f£f /5<f as* 



figure 16.- Shear flow in the case of pure 

bending couples at spar chords. 
Total tending moment M B • 2714 mkg [235,500 
in. -lb.] heavy curve; shear flow determined 
from strain measurement*. Light horizontal 
lines: calculated ▼allies 
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figure 18.- Compar- 
ison of 
calculated and meas- 
ured longitudinal 
stresses in the caee 
of bending by trans- 
verse forces* Load- 
ing case a; Trans- 
verse load 1700 
kg [3760 lb.) at the 
free ends of the 
webs. Loading case 
1 b: Transverse load 
q*340 kg [3M50 lb.] 
on trans veree frames 
b.c f d,e»and f. 



edcbo 'jfedcba 
According to elementary bending theory 
Statically indeterminate calculation 
Averages of stress measurements referred 
to flange centroide 
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Figs. 17,19,20,22 



/ W J&7 W 




I i 




"a 




Figure 19. ~ Decrease of the 

equilibrium force 
groups Tfx introduced by the 
distributed_forces of the 
redundant© X^. 



Figure 17 # - Shear flow in the case of 
bending by a couple at the 
middle of the wing (moment Mb * 2714 mkg 
[235,500 in.~lb.] heavy curve: shear 
flow determined from the strain measure- 
ments. Light horisontal lines: calculat- 
ed values). 



Single - SUSS* + ^llibriun 

«~t.-jj22^+ group 

jfe^J -<L~j3 + Lj=J 

Figure 20.- Component • of the couple « introduced 
into the epar chords (total bending 
moment lf B * 2714 mkg [235,500 in. -la.]). 



, Xlementftrr 
Single _ dlBtribtt ted + Equilibrium 
momenta momeat group 

Figure 22.- Components of the couple introduced 

at the middle of the wing (moment 
M B = 2714 mkg [235,500 in.-lb.]). 
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^^"1 I — - According to elementary bending 
theory 

~) Statically indeterminate 
calculation 
Averages of stress measurements 
referred to flange centroids 



Figs. 21,23 
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(<t) Calculated by means of orthogonal characteristic force groups for 

sheets with equal stiffeners. 
(0) Calculated by means of exact orthogonal characteristic force 
groups of the model. 
Figure 21.- Comparison of the calculated and measured longitudinal stresses 

in the case of couples introduced in the spar chords. 
Loading case c: total moment M B = 2714 mkg [235,500 in. -lb.] 
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1~ "»« lucuBurua longiviiainai stresses 

T^.,. „ _ ™ ^ e ca8e of * couple introduced at the middle of the wine. 
Loading case d: Moment Ife - 2714 mkg [235,500 in.-lb.] 
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